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Abstract 

We study two CP sensitive triple-product asymmetries for neutralino production 
e+ e~ — > ™d the subsequent leptonic two-body decay ^ ^ — > Xi ^) for 

i = e, ^, T. We calculate the asymmetries, cross sections and branching ratios in the 
Minimal Supersymmetric Standard Model with complex parameters /i and Mi . We 
present numerical results for the asymmetries to be expected at a linear electron- 
positron collider in the 500 GeV range. The asymmetries can go up to 25%. We 
estimate the event rates which are necessary to observe the asymmetries. Polarized 
electron and positron beams can significantly enhance the asymmetries and cross 
sections. In addition, we show how the two decay leptons can be distinguished by 
making use of their energy distributions. 
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1 Introduction 



In the Standard Model there is only one physical CP-violating phase in the Cabibbo- 
Kobayashi-Maskawa matrix. The Minimal Supersymmetric Standard Model (MSSM) [1] 
contains several new sources of CP violation if the parameters of the model are complex. In 
the neutralino sector of the MSSM these are the U{1) and SU (2) gaugino mass parameters 
Ml and M2, respectively, and the higgsino mass parameter /i. One of these parameters, 
usually M2, can be made real by redefining the fields. The non- vanishing phases of Mi 
and cause CP-violating effects already at tree level, which could be large and thus 
observable in high energy collider experiments [2]. 

In this paper we study neutralino production (for recent studies with complex param- 
eters and polarized beams see [3, 4, 5]): 

e+ + e- ^ x° + X; (1) 

with longitudinally polarized beams and the subsequent leptonic two-body decay of one 
of the neutralinos 

X- ^ ^+^1, (2) 
and that of the decay slepton 

i ^ X? + 4; i^e,f,,T. (3) 

T-odd observables [3, 6, 7, 8, 9] are a useful tool to study the influence of the CP- 
violating parameters Mi and fi. For the neutralino production (1) and the two-body 
decay chain of the neutralino (2) and (3) we introduce the triple-product 

Ti = (pe- xp^.) -Pi^, (4) 

and define the corresponding T-odd asymmetry 

_ a{Tj > 0) - a{Tj < 0) 

a(T,>0) + a(T,<0)' ^""^ 

where cr is the cross section (33) for reactions (l)-(3). Since under time reversal the triple- 
product changes sign, the asymmetry is a T-odd asymmetry. With the leptonic two-body 
decay of the slepton Eq. (3), we can construct a further T-odd observable which does 
not require the identification of the neutralino momentum. We replace the neutralino 
momentum p^. in Eq. (4) by the lepton momentum pg^ from the slepton decay, which 
defines the triple-product: 

Tii = {Pe- X pe^) -pe^ (6) 

and the T-odd asymmetry 

_ a{Tu > 0) - ajTjj < 0) 
-^'^ a{Tn > 0) + a{Tn < 0) " ^ ^ 
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Due to CPT invariance these T-odd asymmetries are CP-odd if the widths of the ex- 
changed particles and final state interactions are neglected, which is done in this work. 

The T-odd observables in the production of neutralinos at tree level are due to spin 
effects. Only if there are CP-violating phases in the neutralino sector and if two different 
neutralinos are produced, each of the produced neutralinos has a polarization vector with 
a component perpendicular to the production plane [3, 4, 5]. This polarization leads to 
asymmetries in the angular distributions of the decay products, as defined in Eq. (5) and 
(7). In general also spin-spin correlations of the neutralinos are present [5]. However, we 
will not study these in the present work. 

In Section 2 we present the formalism used. In Section 3 we discuss the qualitative 
properties of the asymmetries. We present detailed numerical results in Section 4. Section 
5 contains a short summary and conclusion. 

2 Definitions and formalism 

In this section we give the interaction Lagrangians, the complex couplings and the formu- 
lae for the cross section for neutrahno production (1) and decay (2), (3). For the definition 
of the angles for production and decay, see Fig. 1. 




Figure 1: Definition of angles and momenta of the production and decay process. 
2.1 Lagrangian and couplings 

The interaction Lagrangians for the processes (l)-(3) are (in our notation and conventions 
we follow closely [1, 5]): 

^z^^^ = + Oy^lx?, (8) 
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R?, 



-tlx. 
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-Zj^''[LePL + RePR]£. 



cos 9w 

In the neutralino basis ^,Z, H^, the couphngs are: 



(9) 
(10) 



fii = 



1 



(Tg^ - sin^ 9w)Ni2 + sin 6lTy-/Va 



.cosily 

/^f = -V2e^sin^^[tan^vyiV* -^'a 



O'i- = (^^3A^;3 - Ni,N*^) cos 2/3 - ^ (TV.gAj^ + N^^N*^) sin 2/3, 
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with i, j = 1, . . . ,4. 



■31 



. sin=^ 9 



w- 



(11) 

(12) 
(13) 



where Pl,r = |(1 =F 75)) 9 is the weak couphng constant {g = e/ sin^vK? e > 0), and eg 
and T^e denote the charge and the third component of the weak isospin of the lepton i, 
tan (5 = V2/V1 is the ratio of the vacuum expectation values of the two neutral Higgs fields. 
Nij is the complex unitary 4x4 matrix which diagonalizes the neutral gaugino-higgsino 
mass matrix Y^p, N*^^pN^p^. = m^oSik- 



For the neutralino decay into staus x° 
write for the Lagrangian [11]: 



J^rfxi = 9rkf{aliPR + blPL)Xi + h.c. , A; = 1, 2; i = 1, . . . , 4, 
with 



ffeT, we take stau mixing into account and 

(14) 



Jtj 

"'T3 
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"'TJ 

fR 

Jtj 



with TZl^ the stau mixing matrix defined below and 



R \* 



-Yr{N*^cosP + N*^smp), 



m. 



■/{\/2mw cos/3), 



(15) 
(16) 



(17) 
(18) 



with mw the mass of the W boson and the mass of the r-lepton. The masses and 
couplings of the r-sleptons follow from the hermitian 2 x 2 tl — Tr mixing matrix: 



/ 



— i^L^ ^r) 



M? 



2 I 

tlr I 



TRH 



\^R J 



(19) 
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with 



ML = (ML)* ^mr(A,- I,* te^nP), 
<^f = arg[AT- — /i*tan/3], 



(20) 
(21) 
(22) 

(23) 



At- is the complex trihnear scalar coupling parameter. The f mass eigenstates are (fi, f2) = 
{TL.TRjTV^ with 

f e^'^^cosOf sin^f \ 
7^" = , (24) 

\ — sin 9f e"**^^ cos 6f J 

and 

COS^f = 



-\M? I 



sin = 



ML - < 



The mass eigenvalues are 

= \ + ML) T V(M^.-M4J^ + 4|M^J 

In order to reduce the number of MSSM parameters, we assume 
equations (RGE) for the slepton masses [12]: 

= mo + O.23M2 - m| cos 2(5 sin^ 



(25) 



(26) 

the renormalization group 



= mo + 0.79M| + m| cos W{-\ + sin^ ^w) , 
where is the mass of the Z boson and tuq is the scalar mass parameter. 



(27) 
(28) 



2.2 Cross sections 

In order to calculate the amplitude squared for the complete process of ncutralino pro- 
duction (1) and the two-body decay chain of the neutralino x° (2)-(3), we use the spin 
density matrix formalism of [5, 10]. The amplitude squared can be written as: 

\T\' = 2 $:|A(x°)nA(£)|V(x°)'^Vi(x°)A^A,^2(^~), (29) 

AiA' 

with pp{xf) the spin density production matrix of neutralino Xi , the propagator A(x°) = 
l/[s^^ — m^. + im^Txi], the decay matrix pniiXi) fo^' decay (2), the slepton propagator 
A{£) = l/ls^ — rrij + im^T^], and the factor D2{£) for the slepton decay (3). In Eq. (29), Aj 
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and are the helicities of the neutrahno Xi ^iiid s^. = = pj), the mass and width 

of Xi (^) are denoted by m^^ {mi) and F^. (F^), respectively. The factor 2 is due to the 
summation of the hehcities of the second neutrahno x^, whose decay is not considered. 

The spin basis vectors s^. (a = 1, 2, 3) of the neutrahno Xi fulfiU the orthonormahty 
relations s^. -s^^ = —S"''' and s"^ -p^. = 0. Their explicit form is given in Appendix A. The 
(unnormalized) density matrices can then be expanded in terms of the Pauli matrices: 

ppm'''^ = KKP+j:<ypp^ (30) 

a 

Pdi{x°)kx. = 5a^a.A + E^aja.S^i- (31) 

a 

With our choice of the spin vectors, is the longitudinal polarization of neutrahno 

Xi, is the transverse polarization in the production plane and is the polarization 
perpendicular to the production plane. The analytical formulae for P, Di, D2 and Sp, 
are given in Appendix C. Inserting the density matrices (30) and (31) in Eq. (29) leads 
to: 

|T|2 = 4 |A(x°)r |A(£)|2(PDi + EpEci)L>2. (32) 

The cross section and distributions in the laboratory system are then obtained by in- 
tegrating iTp over the Lorentz invariant phase space element (iLips(s,]9T^ ,p£^,p^j,P£2) 
(B.l): 

o?(7= ^|r|2(iLips(s,p^.,p£i,p^i,p^J, (33) 
where we use the narrow width approximation for the propagators. 

— * — * 

The contributions of the spin correlation terms EpE/^i to the total cross section vanish. 
Their contributions to the energy distributions of the lepton ^l and £2 from decay (2) and 
(3) vanish due to the Majorana properties of the neutralinos [13] if CP is conserved. 
If CP is violated, they vanish to leading order perturbation theory [13]. In our case, 
the contributions can be neglected because they are proportional to the widths of the 
exchanged particles. 

3 T-odd asymmetry 

The asymmetry Ai, defined in Eq. (5), can be written in terms of the angular distribution 
of the decay lepton £1: 

/i° T^^^ cos Oi - L-^ -r^d cos 61 N.-N 
^ ~ Ii T^dcos 01 + ^r^d COS 0i~ N+ + N^' ^ ' 

where COS dj = Fj^ ^^-^^ , ■ -i^^ and thus At is the difTcrcncc of the number of events with 

^ \Pe->^PXi\ \Pll\ 

lepton £1 above (A^+) and below (iV-) the production plane, normalized by the total 
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number of events. In order to measure the asymmetry Ai, the production plane and thus 
the momentum p^. of neutrahno Xi has to be reconstructed. We will discuss in Section 4.4 
how to distinguish the two leptons. 

Analogously the asymmetry An, defined in Eq. (7), can be written in terms of the 
angular distribution of the decay leptons £i and £2'- 

Ii cos en - J-^ -rodeos On 

~ Ii COS On + /o"^ ^r^dcos On ' 

where cos9n = ^^-^^\ ■ r^- Inserting the cross section (33) in the definitions of the 
asymmetries (5) and (7) we obtain: 

_ /Sign[r,,,,]|r|^dLips _ /Sign[r,,,,]Sl,Sl,,dLips 
/irpdLips ~ JPDidLips ' 

with dLips = dLips(s,p^,,P;^.)dLips(s;^,,p^-,p^JdLips(s^-,P;^i,P£2)5(s;^, - m\^?){si - mf), 
sec Eq. (B.l). In the numerator only the spin correlation terms perpendicular to the 
production plane SpS|,^ remain, since only Sp contains the triple-products (4) or (6). 
Thus, the contributions to Ai^n directly stem from Ep. 

In case the neutralino decays into a scalar tau, we take stau mixing into account and 
the asymmetries are reduced due to their qualitative dependence on the x°-ffe-T couplings: 

If 12 _ iLf 12 

M,n oc ' r r ' (37) 

which can be seen from the expressions of Di (C.15) and (C.16). Because the 

asymmetry is only proportional to the absolute values of a^^, 6^^, it is not sensitive to 
a possible phase (pA-r of A^-. In order to be sensitive to 99^1^, one would have to consider 
an asymmetry which involves the transverse polarization of the r. 



4 Numerical results 

In the following numerical analysis we study for = 500 and longitudinally polarized 
beams with P_ = 0.8 and P+ = —0.6, the dependence of the neutralino production cross 
sections cr(e+e~ — >• x°x°), the branching ratios Pi?(x° l€) and the asymmetries Ai 
and An on the parameters \x = \\x\ e^'^'', Mi = \Mi\ e'*'^! and M2 for tan (3 = 10. In order 
to reduce the number of parameters, we assume |Mi| = 5/3 tan^ ^vK-^2 and in Eqs. (27) 
and (28) we take mg = 100 GeV for the slepton masses. Since the pair production of 
equal ncutralinos is not CP sensitive, we discuss the lightest unequal pairs x?X2) X1X3 
and X2 xl- 
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4.1 Production of Xi X2 

In Fig. 2a we show the cross section for X1X2 production for = and ifMi = 0.5 tt in the 
|/i|-M2 plane. The cross section reaches values up to 300 fb. For < 250 GeV the right 
selectron exchange dominates so that our choice of polarization P_ = 0.8 and P+ = —0.6 
enhances the cross section by a factor as large as 2.5 compared to the unpolarized case. 
For > 300 GeV the left selectron exchange dominates because of the larger X2 ~ 
coupling. In this region a sign reversal of both polarizations would enhance the cross 
section by a factor between 1 and 20. 

The branching ratio for the neutralino two-body decay into right selectrons and smuons 
BR(x2 ^di) (summed over both signs of charge) is shown in Fig. 2b. The branching 
ratio reaches values up to 64% and decreases with increasing |//| when the two-body decays 
into the lightest neutral Higgs boson and/or the Z boson are kinematically allowed. 
The channels into the W boson do not open. With our choice for the slepton masses, 
Eqs. (27) and (28), the decay into left selectrons and smuons can be neglected because 
these channels are either not open or the branching fraction is smaller than 1%. As we 
assume that the squarks and the other Higgs bosons are heavy, the decay into the stau is a 
competing channel, which is discussed below. In our scenario this decay mode dominates 
for Af2 < 200 GcV, see Fig. 4a. The resulting cross section a{e^e~ — >■ x?X2) ^ BR(x2 
idi) X BR(4 ^ x%) with BR(£r ^ x%) = 1 is shown in Fig. 2c. 

Fig. 2d shows the |/i|-M2 dependence of the asymmetry An for 99^1 ~ 0.5 tt and 
(fin — 0. In the region |//| < 250 GeV, where the right selectron exchange dominates, 
the asymmetry reaches 9.5% for our choice of beam polarization. This enhances the 
asymmetry up to a factor of 2 compared to the case of unpolarized beams. With increasing 
1/^1 the asymmetry decreases and finally changes sign. This is due to the increasing 
contributions of the left selectron exchange which contributes to the asymmetry with 
opposite sign and dominates for |//| > 300 GeV. In this region the asymmetry could be 
enhanced up to a factor 2 by reversing the sign of both beam polarizations. 

The sensitivity of the cross section a and the asymmetry An on the CP phases can 
be seen by contour plots in the (p/x-ifiMi plane, for = 240 GeV and M2 = 400 GeV 
(Fig. 3). In our scenario the variation of the cross section (Fig. 3a) is more than 100%. In 
addition to the CP sensitive observables, the cross section may also serve to determine the 
values of the phases. Using unpolarized beams, the cross section would be reduced by a 
factor 0.4. The asymmetry An (Fig. 3b) varies between -8.9% and 8.9%. It is remarkable 
that these maximal values are not necessarily obtained for maximal CP phases. In our 
scenario the asymmetry is much more sensitive to variations of the phase ipui around 0. 
On the other hand, the asymmetry is rather insensitive to (pn- For unpolarized beams 
this asymmetry would be reduced roughly by a factor 0.33. 

The relative statistical error of each asymmetry A can be calculated to 6 A = AA/A = 
S/ (AVN), with S standard deviations, assuming a Gaussian distribution of the asymme- 
try A. Here, N — La is the number of events with L the total integrated luminosity and 
a the total cross section. Assuming bA 1, it follows S ^V^- For example, in order 
to measure an asymmetry of 5% with S=2 (confidence level of 95%), one would need at 
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least 1.5 X 10^ events. This corresponds to a total cross section for reactions (l)-(3) of 3.1 
fb with £ = 500 fb~^. We show the contour hnes of 5" = 3 and 5 for An in Fig. 3c with 
jC — 500 fb~^. In the gray shaded area S <3. 

In Fig. 3d we also show the asymmetry Ai which is a factor 2.9 larger than An, because 
in All the CP-violating effect from the production is washed out by the kinematics of the 
slepton decay. However, for a measurement of Ai the reconstruction of the momentum 
is necessary. The asymmetry Ai shows a similar dependence on the phases as An because 
both are due to the non vanishing neutralino polarization perpendicular to the production 
plane. 

It is interesting to note that, due to the weak dependence on the asymmetries can 
be sizable for pa 0. Small values for ^p^ are suggested by constraints on electron and 
neutron electric dipole moments (EDMs) [14] for a typical SUSY scale of the order of a 
few 100 GeV (for a review see, e.g., [15]). 

Next we want to comment on the neutralino decay into the scalar tau and discuss the 
main differences from the decay into the selectron and smuon. In some regions of the 
parameter space, the decay of the neutralino into the lightest stau fi may dominate over 
that into the right selectron and smuon, and may even be the only decay channel. In 
Fig. 4a we show the branching ratio BR{x2 ~^ ^i''') the |//|-M2 plane for A-^ — —250 
GeV, (fiMi = 0.5 TT and (pfj, — 0. For M2 < 200 GcV the branching ratio BR(x2 fir) 
is larger than 80%. However, due to the mixing in the stau sector the asymmetry An, 
Fig. 4b, is reduced compared to that in the selectron and smuon channels, see Fig. 2d. 
The reason is the suppression factor ( | a^j ^ — 1 P)/( | a^j ^ + 1 P) , Eq. (37) , which may be 
small or even be zero. This may lead to a reduced or vanishing asymmetry, respectively, 
even in the case of non zero CP phases. 

4.2 Production of X1X3 

We show in Fig. 5a and b contour plots of the cross section a{e^e~ — > X1X3) ^ BR(x3 — > 
inii) X BR{iR — > Xi^2) with BR(£r — > x?-^2) = 1 and of the asymmetry An, respectively. 
The cross section with polarized beams reaches more than 100 fb, which is up to a factor 
2.5 larger than for unpolarized beams. The asymmetry An, shown in Fig. 5b, reaches 
-9.5%. For unpolarized beams this value would be reduced by a factor 0.75. For our choice 
of parameters the cross section and the asymmetry for XiX^ production and decay show 
a similar dependence on M2 and as for X1X2 production, however, the kinematically 
allowed regions are different. We also studied the dependence of An. For = 
0.57r(0.l7r) and ipMi = 0, the maximal values of An in the M2-\n\ plane are \An\ < 
3%(1%). 

4.3 Production of X2X3 

The production of the neutralino pair e'^e~ —>■ X2X3 could make it easier to reconstruct 
the production plane because both neutralinos decay. This allows one to determine also 
asymmetry Ai, which is a factor 2-3 larger than An. We discuss the decay of the heavier 
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neutralino Xs: which has a larger kinematically allowed region in the |/i|-M2 plane than 
that of Fig- 6 we display the production cross section a{e'^e~ — > X2X3) which 

reaches 100 fb. The cross section a{e+e- X2X3) x BR{xl IrIi) x BR{iR x%) 
with 'BK{iR — ^ Xi^2) = 1 is shown in Fig. 6b. The asymmetry An is shown in Fig. 6d. 
As to the ip^ dependence of Ai, we found that for (^^ = 0.57r(0.l7r) and (^Mi 0, \Ai\ can 
reach 25% (2%) in the M2-\ii\ plane. 



4.4 Energy distributions of the leptons 

In order to measure the asymmetries Ai (5) and An (7), the two leptons li and £2 from the 
neutralino (2) and slepton decay (3) have to be distinguished. We therefore calculate the 
energy distributions of the leptons from the first and second decay vertex in the laboratory 
system (i.e. the cms of the incoming e+ and e" beams). One can distinguish between the 
two leptons event by event, if their energy distributions do not overlap. If their energy 
distributions do overlap, only those leptons can be distinguished, whose energies are not 
both in the overlapping region. 

The energy distribution of lepton in the laboratory system has the form of a box 
with the endpoints: 



mr 



£1 ,min,max 



2(^x.±9)' 



(38) 



with q the neutralino momentum. The energy distribution of the second lepton £2 is 
obtained by integrating over the energy E^ of the propagating slepton: 



1 da 
a dEi^ 

with: 



m| ml 



q[ml — mj][mj — 



X 



A 
a 

b 



<Ee, < 



a 
b 

B 



(39) 



A,B 
a, b 



E. 



i.max.min 



2 

mi 



m 



Xi 



2mi 



E,. 



I, max 



711; 



171; 



m 



XI 



2mi 



ExArrij + m|) ± (ml - m])^^ 
2mJ 

Xi 



(40) 
(41) 

(42) 



We show in Figs. 7a - c an example of the energy distributions of lepton li (dashed 
line), and lepton C.2 (solid line), i = e,fi, for e~^e~ X1X2 subsequent decays 

X2 ik and I x%, for tan/? = 10, M2 = 300 GeV, v?^ = and (fM, = O.Svr 
for 1/^1 = 200, 300 and 500 GeV, respectively. The parameters are chosen such that 
the slepton mass m^^ = 180 GeV is constant, the LSP mass m^^ — 140, 145, 150 GeV 
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is almost constant whereas the neutrahno mass m^^ = 185; 240, 300 GeV is increasing. 
The mass difference between in and Xi decreases (Am = 40,35,30 GeV), whereas the 
mass difference between X2 '^R increases (Am = 5,60, 120 GeV). The endpoints of 
the energy distributions of the decay leptons depend on these mass differences. Thus, in 
Fig. 7a, the second lepton is more energetic than the first lepton. The energy distributions 
do not overlap and thus the two leptons can be distinguished by measuring their energies. 
This also holds for Fig. 7c, where the first lepton is more energetic than the second one. 
In Fig. 7b the two distributions overlap because the mass differences between Xi: and 
X2 are similar. One has to apply cuts in order to distinguish the two leptons, which reduce 
the number of events. 

A potentially large background may be due to slepton production e^e~ — > — > 
i'^i~XiXi- However, these reactions would lead generally to "two-sided events", whereas 
the events from e+e~ — > XiXi ~^ •^'^"•^"XiXi are "one-sided events". Moreover, the back- 
ground reaction is CP-even and will not give rise to a CP asymmetry, because the sleptons 
are scalars and their decay is a two-body one. 

5 Summary and conclusion 

We have considered two CP sensitive triple-product asymmetries in neutralino production 
e"^e~ — > XiX'j and the subsequent leptonic two-body decay chain of one neutralino x° — > 
££, £ — > Xi'^ for = e, T. The CP sensitive contributions to the asymmetries are present 
already at tree level and are due to spin effects in the production process of two different 
neutralinos. The asymmetries are induced only if CP-violating phases of the gaugino and 
higgsino mass parameters Mi and/or /i are present in the neutralino sector of the MSSM. 

In a numerical study for e'^e~ X1X2 and neutralino decay into a right slepton 
X2 ~^ ^R^ ^ave shown that the asymmetries can be as large as 25%. They can be 
sizeable even for a small phase of /x, which is suggested by the experimental limits on 
EDMs. The asymmetries are simlar for the processes e^e~ X1X3 and e+e~ — » X2X3- 
Depending on the MSSM scenario, our proposed asymmetries should be accessible in 
future electron-positron linear collider experiments in the 500 GeV range. Longitudinally 
polarized electron and positron beams can considerably enhance both asymmetries and 
production cross sections. 
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Appendix 



A Momentum and spin vectors 

We choose a coordinate frame in the center of mass system such that the momentum 
of neutraUno x° points in the ;i;-direction. The scattering angle is d^ipe-,Pxj) the 
azimuth can be chosen to zero. The momenta are given by: 

Pe- = — sin^,0,cos6'), Pe+ = -£'6(1, sin ^, 0, — cos^), (A.l) 
Px. = (^x.> 0, 0, -q), p^. = (E^.A 0, q), (A.2) 

with the beam energy Eb = \/i/2 and 



^ ^ ^1, £; ^ ^, q= ^ ' (A.3) 

^' 2^5 ' ' ^ 2^5 



with m^. , m^^. the masses of the neutrahnos and X{x, y, z) — x'^+y'^ + z^ — 2xy — 2xz — 2yz. 
The three spin basis vectors of x° are chosen by: 

4 = (0,-1,0,0), 4 = (0,0,1,0), 4 = ^(g,0,0,-Ej. (A.4) 

Together with the unit momentum vector Vxil'^Xi , the spin basis vectors form an 
orthonormal set. The momenta and energies of the leptons are: 

= (E^j, sin6'i COS01, Ip'^J sin6'i sin0i, -|p^J cos6'i), (A. 5) 

= {Eii^ , - |p^2 1 sin ^2 cos 02, \Vh I sin sin 02, - 1 cos 612) , (A. 6) 

2(^xi -^cos^i) 2(E^-- Ipx^ J cos ^1^2) 

with 6*1 = L{pi^,p^X 02 = l{pe^,p^i) and the decay angles OD^^iPhP^i) ^nd 9Dil{Pxi,Pi) 
(see Fig. 1): 

cos^iij = cos^Dj cos^2 — sin^£)^ sin^2cos(02 — 0i), cos^j)^ = — ^t^"- (^-8) 

IPxi I \Pxi - Pil I 

B Phase space 

The Lorentz invariant phase space element for the neutralino production (1) and the decay 
chain (2)-(3) can be decomposed into the two-body phase space elements: 

dUps{s,Py^.,Pe„p^,,Pi,) = 

dUps{s,p^.,p^.) ds^. dUps{s^.,pi,ptJ ds^ dU^s{si,p^^,pt^), (B.l) 
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dUps{si,p^,,Pe2) 



sin 6* d6, 



2(27r)2 ml. - mj 



1 



dQi, 



2(27r)^ mj — m 



dn.2, 



XI 



with s^. — p^^, = pj and dJlj = sin^j d9i d(f)i. 



(B.2) 
(B.3) 
(B.4) 



C Neutralino production and decay matrices 

In this Section we give the analytical expressions for P, Sp and Di,D2, "^di in the center 
of mass system. Expressions for Ep^ can be found in [5] . 

C.l Neutralino production 

The analytic expression P of Eq. (31) is independent of the neutralino polarization. It 
can be decomposed into contributions from the different production channels [5] : 



P = P(ZZ) + P(ZeR) + P(ZeL) + Pienen) + P(eLei), 
with 



(C.l) 



P{ZZ) 



P{ZeR) = 



9^ 



COS^ U]v 



A^iZ)m - PtPtm + Ri) - {Pi - Pi){Li - Ri)]Ei 



{\0:f\'{E^^E^^+q'cos'e) - [{ReOl^r - (/mOy)>,,m,J, 



(C.2) 
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cos^ 6w 



Reil + P!)(1 - Pl)ElRe{A'{Z) 



- {^'\en)f^*fp:^* + A«*(g^)/«/5*0:;^)m,,m 
-(A**(e~«)/5750y - A-*{^.n)ifj^*0:^*)2E,qcose 
+iA'*{en)f^*hW + A"*(ei^)/e1/et0y*)(^x.^x. + cos^ 9)]}, 
l(l + P^_)(l-Pl)E!{\f^f\l 



(C.3) 



R\2 

ej\ 



{\A\eR)\'+\A^{iR)\'){E^^E^^+q'cos'e) - {\A\eR)\' - \A\iR)\')2E,qcose 
- i?e{(/,t)2(/5)2A«(e^)A^*(e«)}2m,,m,J. (C.4) 

To obtain the quantities P{Zei), P{eLei) one has to exchange in Eqs. (C.2) - (C.4) 



A'ien)^ A\iL), A"(gK) ^ A"(gi), P! ^ P^, P^ ^ P^ 



Rp 



R 



o. 



fR ^ f 

J ei J eii J t 



(C.5) 
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The propagators are defined as follows: 



s — niz + tnizL z t — m^^^^ u — rn^^^ 

where m and V denote the mass and width of the exchanged particle, respectively, and 
s = {pe- +Pe+)^, t = (Pe- ~PxjYy u = {pg- —PxiT- The longitudinal beam polarization 
of e~{e^) are denoted by Pl{Pl), respectively. Generally the contributions from the 
exchange of cr (e^) selectron exchange is enhanced and that of cl (e^) is suppressed for 
> 0, P| < (Pf < 0, P3 > 0). 

C.2 Neutralino polarization 

The analytic expressions for the coefficient Sp in Eq. (30) which describes the transver- 
sal polarization of neutralino x° perpendicular to the production plane decomposes into 
contributions from the different production channels [5]: 

^piXi) = ^pixl ZZ) + EUxl Zen) + E?,(x°, Zcl) + E|>(x°, e^e^) + E?,(x°, eie,,).(C.7) 
In the center of mass system they read [5] : 

^l{xlZZ) = -4{^^)'\A^{zmi-PlPl){Ll-Rl)-{Pl-Pl){Ll + Rl)] 

^ cos u\Y ' 

X m^.qEl sin QRe{d[f)Im{d[f), (C.8) 

Y?p{x\.Z~en) = -^^Re{l + P'^) {I ^ PDm^.E^q sin 9 

cos U]y 

X Im{A%Z)[f^J,fo:fA^*{en) - f^*fep:f*^'*{en)]}, 

(C.9) 

^Uxlenen) = -^{1 + P'_){1 - Pl)m^^E',qsme 

X /m{(/,t)2(/5)2A"(e«)A**(e«)}. (C.IO) 

To obtain the expressions for T?p{^, Zcl) and Ep(x^,eLeL) one has to apply the ex- 
changes (C.5) in Eq. (C.8) - (C.IO). 

C.3 Neutralino decay matrix 

The neutralino decay matrix is given by Eq. (31): 

a 

The expansion coefficients Di and E^^^ for the decay into the right slepton are (£ = e, /x): 

A - yi/.frK-^f), (c.ii) 
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respectively, with s^. the neutrahno spinvector and the lepton ii momentum vector. 
For the decay into the left selectron or smuon they are: 




For the decay into the stau (A; = 1, 2) they read: 




The decay of the right {R) or left (L) slepton into lepton and xl is given by {i — e, /i): 
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Figure 2: Contour plots for 2a: cx{e+e- x?^"). 2b: BR(x^ i = e,fi, 2c: 

^(e+e- ^ x?X2) X BR(x^ ^ ^Vi) x BR(£r ^ x?^2) with BR(4 ^ X?^2) = 1, 2d: the 
asymmetry An, in the \fi\-M2 plane for 92^/1 = O.Svr, = 0, taking tan/5 = 10, uiq = 100 
GeV, = -250 GeV, = 500 GeV, P_ = 0.8 and P+ = -0.6. The area A (B) is 
kinematically forbidden by m^o + m^o > ^/s (?^^^ > ^^o)- The gray area is excluded by 
m-± < 104 GeV. 

Al 
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a(e+e-^X?X?4^2) in fb 
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Figure 3: Contour plots for 3a: (7(e+e- XiXl) x BR(x2 ^ ^di) x BR(£r x%) 
with BK{iii ^ xV'^) — 1) 3b: the asymmetry An, 3c: the standard deviation S, 3d: 
the asymmetry Ai, in the ip^-ipM^ plane for M2 = 400 GeV and = 240 GeV, taking 
tan/5 = 10, mo = 100 GeV, = -250 GeV, ^ = 500 GeV, P_ = 0.8 and P+ = -0.6. 
For ipMi,^ij. — we have m^^ = 221 GeV, m^o — 178 GeV and m^o — 243 GeV. 
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Figure 4: Contour plots for 4a: BR(x2 ^i^) and 4b: the asymmetry An, in the \fi\-M2 
plane for ipMi = O.Svr, = and Ar = —250 GeV, taking tan/3 = 10, mo = 100 GeV, 
^/s = 500 GeV, P_ = 0.8 and P_ = —0.6. The area A (B) is kinematically forbidden by 
m^o + m^o > ^/s {nif^ > m^o). The gray area is excluded by m-± < 104 GeV. 
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Fig. 5a /GeV Fig. 5b \n\ /GeV 



Figure 5: Contour plots for 5a: a{e+e- x^xl) x BR(x3 ^ x BR(4, X%) 

with BR(£ij — * Xi^2) = 1 and i = e,fi, 5b: the asymmetry An, in the |/i|-M2 plane for 
cpj^^^ = 0.57r, (pf, = 0, taking tan/3 = 10, mo = 100 GeV, = -250 GeV, = 500 GeV, 
P- = 0.8 and P+ = —0.6. The area A (B) is kinematically forbidden by m^o + m^o > ^/s 
[mi > m^o). The gray area is excluded by m-± < 104 GeV. 
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Figure 6: Contour plots for 6a: cr(e+e~ — * X2X3)) 6b: a{e^e^ X2X3) ^ BR(x3 
Inh) X BR(£/j — > Xi^2) for £ = e,/x, and BRilu — * Xi^2) = 1, 6c: the asymmetry Ai, 
6d: the asymmetry An, in the |yu|-M2 plane for ipM^ = O.Svr, ip^ = 0, taking tan/3 = 10, 
mo = 100 GeV, = -250 GeV, = 500 GeV, P_ = 0.8 and P+ = -0.6. The area 
A (B) is kinematically forbidden by m^o + m^o > i^iji > "^xjj)' "^^^ S^ay area is 
excluded by m-± < 104 GeV. 

Al 
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Fig. 7c E /GcV 

Figure 7: Possible types of energy distributions in the laboratory system for ii (dashed 
line) and ^2 (sohd line) for e^e'^ X1X2 ^^"^ subsequent decays X2 ^nh and 
h X%: for M2 = 300 GeV, m^^ = 180 GeV, tan/3 = 10 and {\n\,m^„m^,}/GeV = 
{200, 140, 185}, {300, 145, 240}, {500, 150, 300} in a, b, c respectively. 
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